World citation and collaboration networks: uncovering the role of geography in science 
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Modern information and communication technologies, especially the Internet, have diminished 
the role of spatial distances and territorial boundaries. This has enabled scientists for closer collab- 
oration and internationalization. Nevertheless, geography remains an important factor affecting the 
dynamics of science. Hence, assessing the influence of spatial proximity between scientists is crucial 
to promote efficient collaboration strategies and, ultimately, to improve the quality of science. Here 
we present a systematic analysis of citation and collaboration streams between cities and countries. 
By assigning papers to the geographic locations of their authors' affiliations, we construct weighted 
networks of citations and collaborations. The citation flows as well as the collaboration strengths 
between cities decrease with the distance between them and follow gravity laws with exponents close 
to 1. Moreover, for a given number of authors, the diversity of affiliations increases the number of 
citations, especially when many countries are represented. In addition, the total research impact of 
a country grows linearly with the amount of national funding for research & development. How- 
ever, the average impact reveals a peculiar threshold effect: the scientific output of a country may 
reach an impact larger than the world average only if the country invests more than 120,000 US 
$ per researcher annually. Our results reveal the overall structure of scientific research by showing 
the correlation between collaboration, citation, geography and funding, and could provide valuable 
inputs in shaping the future science policies. 



I. INTRODUCTION 

The strength of most interactions in nature typically 
decreases with the distance between objects or con- 
stituents. The most famous example is Newton's gravi- 
tational force, which is known to decay with the square 
of the distance between the masses. This principle holds 
also outside the realm of physical processes. Recent stud- 
ies on mobile phone communication networks [TJ [2] and 
blogs [3 have revealed that the probability for a social 
tie to occur between agents decays with a power of their 
distance. 

Likewise, scientific interactions are likely to take place 
between scholars localized in the same or nearby areas. 
Scientists tend to cluster in space, since the elaboration 
and progress of a project requires frequent discussions 
between collaborators that is hardly possible if they live 
far apart. Factors based on cultural, linguistic and insti- 
tutional differences cause additional obstacles to long- 
distance cooperation [4|. Further, research funding is 
mostly allocated at the national level [5 , thus favoring 
regional over international collaborations. 

Nowadays, the Internet and the greater affordability 
of international transportation have enormously reduced 
distances between people, overcoming both geographic 
and cultural barriers [6-8 . This in turn has made sci- 
entific collaborations between distant scholars far easier 
than before [9-12 . Multi-authored collaborations serve 
as big opportunity for science [13 j, as one can integrate 
a wide range of competence and skill, to attack diffi- 
cult problems, with an enhanced chance of success. In- 
deed, the last decades have witnessed the formation of 
larger and larger research teams [HI [15]. In particular, 



mult i- university collaborations have been growing at a 
fast pace and are more likely to lead to high impact pub- 
lications [16], especially if they involve different coun- 
tries [HI HE]. 

Geographic proximity is also likely to affect the pro- 
cess of giving and receiving credits for someone's work, 
expressed by paper citations. For most papers one ex- 
pects to find a decaying probability of citation with dis- 
tance, as new findings are typically more visible in the 
area where the authors operate. This is confirmed by a 
recent study [19] . In addition, collaboration patterns are 
likely to influence and be influenced by citations. While 
collaborating, scholars become more familiar with the 
scientific output of their co-authors, which then has a 
higher chance to be cited in the future. In turn, scholars 
citing frequently each other's work have strongly over- 
lapping research interests, and are more likely to become 
co-authors sooner or later. Therefore citations and col- 
laborations between distinct locations are likely to be 
correlated. 

Knowing how scientific interactions vary with distance 
is also valuable for practical reasons. To scholars, it 
might suggest how to choose collaborators in order to 
optimize the impact and visibility of their research. To 
institutions and governments, it might advice suitable al- 
locations of funds for regional and international projects, 
in order to improve the scientific outcome for a given 
amount of resources. It is then not surprising that spatial 
scientometrics has acquired a prominent role during the 
last few years. There are a number of studies carried out 
exploiting the enhanced availability of citation data [20J . 
Yet there are other factors, namely funding, that also 
plays a crucial role in the development of a research 
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Figure 1. Properties of the world citation network. (A) Citation map of the world where the area of each country is scaled and 
deformed according to the number of citations received, which is also represented by the color of each country. (B) Citation 
distribution of papers of top 20 countries. If a paper is written by authors from multiple countries, the paper contributes to 
each country. (C) When the distributions in (B) are normalized by the average number of citations of each country, they fall 
on top of each other. (D) Probability distribution function of the number of citations received by each city. (E) Cumulative 
distribution function of the link weights Wij (excluding self-links) and self-links wu in the citation network of cities. (F) Node 
in-strength against its in-degree for the city citation network. (G) Link weight against the product of the strengths of the 
connected nodes in the city citation network. For each plot we show the corresponding best-fit lines and power law exponents. 



project, as it not only contribute towards the direct and 
overhead costs of the research but also facilitates the co- 
operation and collaboration among researchers working 
in different locations and different fields [2T] . Since both 
public and industrial resources are used to fund academic 
research, it is also natural to question the result and im- 
pact obtained with these resources [22j [23] . 



We have performed the first comprehensive study of 
citation and collaborative interactions between different 
geographic locations. We used one of the world's largest 
citation databases to derive the citation and the collab- 
oration network, i.e. weighted networks where nodes are 
cities and links are citations and collaborations between 
the corresponding cities (see Methods). The analysis of 
these networks [24-27 discloses the existence of gravity 
laws as well as non-trivial correlation between collabo- 
rations and citations. We also present a global study of 
the relationship between funding and scientific outcome 
of different countries. 



II. RESULTS 

The research contribution of each country in terms of 
the (normalized) number of citations received iVcite is il- 
lustrated in the world map of Fig. [TJ\. Colored maps can 
be misleading as the value assigned to a large area gives 
an impression of a much greater impact of that color 
in the visualization. We thus created a cartogram, in 
which the geographic regions are deformed and rescaled 
in proportion to their relative research contribution [28J. 
The citation strengths of countries span over seven or- 
ders of magnitude. North America and Europe receive 
42.3% and 35.3% of world's citations, respectively. In 
contrast, the contribution by Asia amounts to only 17.7% 
of world's citations while the total contribution of Africa, 
South America and Oceania is lower than 5%. In this 
ranking the United States is the leading country fol- 
lowed by the United Kingdom, Germany, Japan, and 
China. The corresponding world map in terms of coun- 
tries' number of (normalized) publications is shown in 
the Appendix, Fig. |B.l^ . 

In Fig. [Tj3 we plot the probability distribution of the 
number of citations of papers in the largest 20 countries. 
A paper is associated to a country if at least one of its 
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affiliations is from that country. All these distributions 
are broad and vary over four orders of magnitude. When 
each distribution is rescaled by the average number of 
citations of papers of the respective country, all curves 
nicely collapse (Fig. [lp). This result suggests that the 
functional form of the citation distribution is the same 
in each country and that the difference between countries 
can be effectively summarized by the average number of 
citations. This type of universality holds at t he level 
of scientific disciplines as well [29 . In Fig. B.2A. of the 
Appendix shows the world map in terms of the average 
number of citations per paper of each country. 

Next we consider the contribution at the level of cities. 
In Fig. [Tp we plot the probability distribution of the 
cities' citations. The distribution is broad, spanning over 
five orders of magnitude, and it follows a power law decay 
with exponent 1.46 ± 0.03. This suggests a relationship 
with the population of the city, as the city size distribu- 
tion obeys the Zipf law [30j [31], i.e. decays as a power 
law (with exponent 2). 

We now consider the weighted citation network be- 
tween cities, where the nodes are the cities that are con- 
nected by weighted and directed links, indicating pub- 
lications of one city citing publications of the others. 
The network has 18,199 nodes and 9,494,021 links in- 
cluding 14,447 self-links (i.e., citations within the same 
city). In Fig. [Tj3 we plot the cumulative distribution 
of the weights of self-links and links between differ- 
ent nodes. Both these distributions are broad; how- 
ever, the weights of self-links are more heterogeneous, 
revealing a bias towards self-citations. Next we calculate 
the number of incoming links, i.e., the in-degree kf 1 of 
each node i and its in-strength, sf 1 = ^jwQ te , which 
equals the number Np lte of (normalized) citations re- 
ceived. By plotting the in-degree against the in-strength, 
we find that there is a power law scaling behavior with 
(s m )(fc m ) oc (k in ) a (Fig.^E). However, there are two dis- 
tinct scaling regimes: for nodes with small kf 1 (< 200) 
the exponent is a = 0.91 ± 03 (regression coefficient 
R = 0.95), while for large kf 1 (> 200) the exponent is 
a = 2.20 ± 0.08 (R = 2.01). The super-linear behav- 
ior suggests that stronger links are more frequently con- 
nected to high in-degree nodes. The out-strength of the 
nodes follows a similar relationship with the out-degree of 
the nodes (see Appendix, Fig. |B.lp ). Finally, we plot the 
weights of the links wff te against the product of the node 



strength s° ut s in . The product s° 



gives the weight of 



a link that is expected to occur by chance between i and 
j if all the papers would be citing each other at random. 
Even in this case there are two distinct scaling regions, 
wg- ite ex (s° ut sf) a , where a = 0.13 ± 0.01 (R = 0.19) if 
the product is less than 2 x 10 T , while for larger values 
of the product a = 0.99 ±0.01 (R = 1.07). This suggests 
that the observed citation is as expected between high 
strength nodes, while it is much lower in case of cities 
with low strength. 

Let us now consider the collaboration network at the 
city level, where the nodes are cities and weighted undi- 
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Figure 2. Properties of the world collaboration network. (A) 
Cumulative probability distribution of the link weights in the 
collaboration network of cities. Self-links are shown sepa- 
rately. (B) Fraction of internal collaboration, indicated by 
the ratio of the weight w^° l of the self-link and strength Si of 
a node, against Si. (C) Strength of a node against its degree. 
The straight line indicates a power law behavior with expo- 
nent 1.66 =b 0.04. In these plots we use the same colorbar as 
in Fig. 1. 



rected links indicate the presence and frequency of col- 
laborations between scholars of different cities. There are 
18,199 nodes in the network and 1,256,718 undirected 
links including 14,954 self-links. The weight of the self- 
links indicates the amount of internal collaboration. The 
degree of a node i indicates the number of other cities 
with which i collaborates and its strength is indicative 
of, but not coincident with, the number of papers writ- 
ten by scholars of institutions in that city. 

In Fig. |2^ we plot the cumulative probability distri- 
bution of link weights. As for citations, the weights of 
self-links are more broadly distributed than the weights 
of the links between different cities, showing that schol- 
ars of a city collaborate more frequently with each other 
than with colleagues from any other city. The distribu- 
tions of collaboration and citation streams between cities 
differ from their analogues in mobile phone communica- 
tions and world trade, that show log- normal distribu- 
tions [21 [32]. Next, we consider the fraction of internal 
collaboration by calculating the ratio of the weight of 
the self-link to the strength of the node. By plotting 
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j S{ against the strength of the node s^, we see that 
the ratio increases with indicating that as the city size 
increases most of its collaborations take place within the 
city (Fig. [2^). However, for small cities most of their pa- 
pers are written with external collaborators. The node 
degree scales with its strength as (s){k) oc k a , where 
a = 1.66 ± 0.04 (R = 1.65) (Fig. |2p). This super-linear 
scaling suggests that higher degree nodes are more fre- 
quently connected by stronger links. 

Let us explore the relationship between the citation 
and the collaboration networks at both the country and 
the city level. At the country level the collaboration net- 
work comprises 226 nodes and 10,308 undirected links, 
including 219 self-links. In the citation network there 
are also 226 nodes but 28,869 directed links, including 
215 self-links. In Fig. [3J we plot the weight of links of 
the collaboration network, wfj° l against the weight of the 
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Figure 3. Correlation between the world citation and collab- 
oration networks. Weight of the links in the citation network 
against the corresponding links in the collaboration network 
at the (A) country level and (B) city level network. Power law 
scaling is shown by solid lines with exponents 1.04 ±0.01 and 
0.82 ±0.02, respectively. Density plot of the number of cita- 
tions of a publication against the number of (C) co-authors, 
(D) countries (E) cities in the affiliation. The circles indicate 
the average trend. 



same links in the citation network, w^ lte + 

where a 



Cite 
ji ' 

scaling w^° l oc (wg ite ± w^ te ) a where a = 1.04 ± 0.01 
(R = 1.04) for countries (Fig. [3^), and a = 0.82 ± 0.02 
(R = 1.05) for cities (Fig. [3^>), i.e. the increase in col- 
laboration is linearly related to the amount of citations 
exchanged between the two countries/cities. 

It has been previously shown that papers published 
by teams often get more citations than single author pa- 
pers p~5| [16] . Our results also show that the average num- 
ber of cites of a publication increases with the number 
of co-authors of that publication (Fig.j3p). Furthermore, 
the average number of citations of a publication increases 
with the number of affiliated countries and cities of its au- 
thors (Fig.|3p and E). In Table[l) we consider papers with 
a given number of authors and categorize them accord- 
ing to whether all the affiliations listed in the paper are 
from a single city, from multiple cities in a single coun- 
try or from different countries. For an equal number of 
authors, publications having multiple international affili- 
ations get a statistically significant increment (p < 10 -4 ) 
in the number of citations with respect to publications 
with only domestic affiliations. Thus, crossing territorial 
boundaries also pays off in terms of scientific impact. In 
contrast, multiple domestic affiliation do not positively 
effect the number of citations when the number of au- 
thors in a publication is less than 6. 



Table I. Dependence of citations on collaboration. We cate- 
gorize each paper by the number of authors and their affili- 
ations. For each of these groups we indicate the fraction of 
papers that are in the group and the mean number of cita- 
tions. The error represents the standard error of the mean, 
calculated using bootstrap sampling with repetition. 
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City 


City 


Countries 


1 


13.03 


4.25 ± 0.02 


4.95 ± 0.12 


5.24 ± 0.11 


2 


19.01 


6.80 ± 0.02 


6.11 ± 0.04 


7.00 ± 0.05 


3 


18.34 


6.92 ± 0.02 


6.38 ± 0.03 


7.30 ± 0.04 


4 


14.95 


7.19 ± 0.02 


7.02 ± 0.03 


8.03 ± 0.04 


5 


11.10 


7.62 ± 0.03 


7.66 ± 0.03 


8.79 ± 0.04 


6 


8.01 


8.13 ± 0.04 


8.52 ± 0.05 


9.77 ± 0.05 


7 


5.20 


8.85 ± 0.05 


9.56 ± 0.07 


10.90 ± 0.07 


8 


3.45 


9.50 ± 0.07 


10.67 ± 0.09 


12.10 ± 0.10 


9 


2.22 


10.23 ± 0.10 


11.52 ± 0.12 


13.17 ± 0.12 


10 


1.53 


10.57 ± 0.12 


12.45 ± 0.14 


14.70 ± 0.15 


>10 


3.17 


13.82 ± 0.17 16.64 ± 0.16 


21.37 ± 0.17 



Next we consider the effect of geographical proximity 
on the citation and collaboration networks by determin- 
ing the geographic location (latitude and longitude) of 
each place in the dataset [33 (see Methods). We found 
that the probability that there is a link between two cities 
in the collaboration network decreases as a power law as 
the distance between the two cities increases (Fig. [4]A.). 
The power law exponent is 0.57 ±0.01. Our results 
are different from those obtained in Ref [34 , where it 
was found that the distribution of distances between co- 
authors decreases exponentially. Such difference might 
be due to the limited dataset used in Ref [34], which in- 
cluded only papers published before 1990, and possibly 
also due to the recent advances in communication and 
transportation technologies. 

Many spatially embedded networks have been observed 
to follow gravity laws [33 , where the flow between two 
locations follows 



PiPi 



(i) 



Here, is the flow between nodes i and j, Pi and Pj 
are the populations of nodes i and j, respectively and 
dij is the Euclidean distance between i and j, the value 
of exponent a being dependent of the system. For the 
collaboration network Eq. [I] becomes 
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In Fig. 4^3, we plot the ratio wf° l /{siSj) against the dis- 
tance afj between all node pairs. We found that as the 
distance increases {wf? l /{siSj)) decreases as a power law 
with the exponent a = 1.16 ± 0.03 (R = —0.97), ex- 
cept at very short distances. As we have seen before, 
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Figure 4. Effect of geographical proximity in the world collab- 
oration and citation networks. The probability of existence of 
a link as a function of the distance between two cities in the 
(A) collaboration network and (B) citation network. Distribu- 
tion of the ratio of the link weight and product of the strengths 
of its endpoints in (C) collaboration network, wf° l /siSj and 



(D) citation network, 



e /s° u s) n against the distance dij 



between the cities. For each distance the average ratio is also 
shown. The solid line indicates a power law behavior with 
exponent a = 1.16 ± 0.03 and 0.77 ± 0.02 respectively. 



collaboration and citation between two places are corre- 
lated. Hence, we also look at the geographical proximity 
in the citation network. We found that the probabil- 
ity that there is a link between two cities in the cita- 
tion network also decreases with distance as a power law 
(Fig. pb). In this case the power law exponent is much 
lower (0.30 ±0.01). The gravity law for the citation net- 
work reads 



Cite S i S 3 



(3) 



In Fig. [ZJl) we plot wfj lte / (s° ut s™) against the distance 
between all the node pairs in the citation network. As for 
the collaboration network we found that (wfj lte / (s^s 1 ^)} 
decreases with distance as a power law with the exponent 
a = 0.77 ± 0.02 (R = —0.35). The above analysis shows 
the existence of an important spatial component in both 
the citation and the collaboration network. It shows that 
both our collaborators and our citations typically come 
from our spatial neighborhood. Further, long distance 
collaborations as well as citations decrease as a power law 
of distance. The difference of the scaling exponents of the 
two networks suggests that two distant places are more 
likely to cite each other than collaborate. Additional 



results are shown in the Appendix, Fig. B.3 
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Figure 5. Relation between research outcome and funding. 
Average number of citations per paper of a country against 
(A) the expenditure in research and development (in millions 
of dollars per year, and purchasing power parity) and (B) the 
number of researchers in that country. The solid line indicates 
power law scaling with exponent 0.99 ± 0.03 and 0.98 ± 0.04, 
respectively. (C) Average number of citations per paper of a 
country against the average spending per researcher. The hor- 
izontal line indicates the average number of citations over all 
papers of all countries, the vertical line indicates the threshold 
of 120,000$ per researcher per year. 



The research performance of each country is gener- 
ally estimated on the basis of the number of publica- 
tions and citations. Although these are straightforward 
measurements of research output, they depend on a wide 
spectrum of resources [35J. For instance, the number of 
researchers and facilities (instruments, laboratories, li- 
braries and other resources) available are typically dif- 
ferent in different countries. A key determinant is the 
funding available for research & development (R&D). To 
quantify the expenses in R&D of a country we consider 
the fraction of gross domestic product (GDP) that is 
spent on R&D. To get rid of economic inequalities in dif- 
ferent countries we consider the R&D spending in terms 
of the purchasing power parity (PPP). In Fig. [5^, we 
plot the number of citations Nca e against the R&D ex- 
penditure and find that it scales linearly with funding. 
Such correlation is not surprising, but the scaling expo- 
nent is non-trivial. It suggests that it is not possible 
to perform or contribute substantially unless there is a 
corresponding amount of funding available for research. 
Moreover, the research contribution in terms of citations 
also scales linearly with the number of researchers in that 
country (Fig. [5^). This result is consistent with the fact 
that the R&D expenditure is correlated with the number 
of researchers. The number of publications of a country 
also shows similar scaling against R&D expenditure and 
number of researchers (Appendix, Fig. B.4). 

Finally as a measure of impact of a country's scientific 
output we consider the average number of citations to the 
publications of that country. In Fig.[5p we plot this num- 
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ber against the average spending per researcher per year 
(R&D expenditure divided by the number of researchers) . 
The latter is not the average salary of researchers in 
that country, as it includes other expenditures such as 
infrastructure, bureaucracy, instruments, etc. This plot 
is much more scattered than the previous plots and does 
not show any definite correlation pattern. Nevertheless, 
the plot reveals a striking feature. If the average spend- 
ing is less than 120,000$ (vertical line in the plot) per re- 
searcher per year we see an increase in the average num- 
ber of citations with spending. However if the average 
spending exceeds this limit, it becomes scattered and in- 
dependent of funding. This figure shows that very rich 
countries like Kuwait and Luxembourg have high fund- 
ing per researcher, still the average number of citations 
per paper is low. Countries like India, Brazil have high 
funding per researcher as well, but low average number of 
cites; this might mean they are investing more on infras- 
tructure. Switzerland, Costa Rica, Panama, Germany, 
Austria, Netherlands, United States have high spending 
per researcher and their average number of citations is 
also high. If we display the number of cites per paper 
averaged over all countries (horizontal line), we see that 
there are no countries in the top left quadrant, i.e. it 
is not possible to do better than the world's average un- 
less there is sufficient spending. Additional measures of a 
country's research performance and corresponding rank- 



ings are reported in the Appendix, Table B.l 



international participation. As a consequence, we expect 
to observe an increasing tendency to form large teams 
with members of many different countries in the future. 

We have also found that in science money acts as a kind 
of threshold motivator: if one does not pay people enough 
they will not be motivated and the outcomes of the re- 
search are poor; if people are paid sufficiently to take the 
issue of money off the table, internationally competitive 
findings are within reach. On the other hand, for con- 
ceptual and creative tasks, paying more than a certain 
threshold does not necessarily increase the output [37] - 
[39] . Further, our analysis reveals that at the country 
level funding has a positive linear impact on the research 
output both in terms of number of publications as well 
as citations. 

In the future we plan to study the role of cities' popu- 
lation, in particular on the distributions of citation and 
collaboration strengths along with their flows. It is well 
known that most characteristics of cities are strongly cor- 
related to the size of their populations [40]. Furthermore, 
an analysis of the evolution of the world citation and 
collaboration networks would show how the spatial di- 
mension of science dynamics has been affected by the 
progress of technology, internationalization and extreme 
events (e.g. wars, economic crises). This way one could 
infer how the scientific landscape has been shaping up 
in the last decades and how is it possible to create more 
efficient partnerships, via dedicated funding programs at 
the national and/or international level, and consequently 
a more productive and successful scholarly world. 



III. DISCUSSION AND CONCLUSIONS 



Our thorough analysis of the world citation and col- 
laboration networks has revealed that the effects of ge- 
ography on the dynamics of science are relevant, despite 
the recent advances in communication and transporta- 
tion. The occurrence of gravity laws for both citation 
and collaboration implies a preference by scientists to 
interact with peers in their geographic areas. However, 
long-distance interactions are not rare, as the interaction 
strength and probability are characterized by power law 
decays. Our work follows similar findings in mobile phone 
communication [TJ [2] , social media [3] and international 
trade [36 , reinforcing the belief that gravity laws hold in 
several different contexts, and that scientific interactions 
are not exceptional from this point of view. 

Citation and collaboration streams between distinct lo- 
cations are strongly correlated, with an approximately 
linear relation. An increase in the number of collabora- 
tions between two cities is then expected to be followed 
by a proportional increase in the flow of citations be- 
tween the cities. This is justified from the fact the peo- 
ple/groups working in similar fields and subject area are 
more likely to cite as well as collaborate with each other, 
and also suggests a natural bias towards self-citation, 
of which we have provided strong quantitative evidence. 
From the point of view of scientific impact, it pays off for 
a team to put together several institutions with a strong 



IV. METHODS 
A. Data description 

We have analyzed all publications (articles, reviews 
and editorial comments) written in English from 2003 
till the end of 2010 included in the database of the 
Institute for Scientific Information (ISI) Web of Sci- 
ence. For each publication we extract the affiliations 
of the authors and the corresponding citations to that 
publication. We parsed the affiliations of all publica- 
tions and have determined the geographic location at 
the city and country level. If there are multiple affil- 
iations listed in a publication, the latter is associated 
with all represented cities and countries. After obtain- 
ing the locations we use the publicly available resources 
(www.wikipedia.org and maps.google.com) to deter- 
mine their coordinates (latitude and longitude). Our 
dataset consists of 8,094,948 publications which have re- 
ceived 62,105,592 citations during the period 2003-2010. 
We were able to extract the geographical information 
from 8,092,314 publications. Affiliations refer to 226 
countries and 37,750 cities. In order to get rid of anoma- 
lies due to any misclassification, we have only consider 
those places that have appeared in at least 5 publica- 
tions during the period 2003-2010. This cutoff led us to 



7 



18,199 cities, producing 99.8% of the total publications 
and receiving 99.9% of total citations. 

Country level information regarding expenditures for 
research and development (R&D) in terms of purchas- 
ing power parity (PPP) and number of researchers 
in R&D are obtained from the World Bank Data 
(databank.worldbank.org) for each year between 2003 
till 2010. By aggregating these yearly datasets we de- 
termine the average of each of the above quantities for 
the period 2003-2010. The data of expenditure for R&D 
is available for 102 countries, the numbers of researchers 
for 89 countries and for 77 countries both informations 
are available. Further details can be found in the Ap- 
pendix [Aj 



B. Network construction 

We have analyzed the data at the country and the city 
level. As the publications and their affiliations form a 
bipartite graph, we construct the collaboration network 
between countries (cities) by projecting it onto the space 
of affiliations. In this collaboration network individual 
countries (cities) act as nodes, and links between them 
indicate that they have appeared in the same publication. 
If a paper is written by authors with n affiliations, we 
put \n x (n — 1) undirected links between each possible 
pair of collaborating countries (cities), with every link 



having weight nx ^_ 1 ) • The total weight between any 
pair of nodes is the sum of all the weights over all the 
publications in the dataset. If there is a single affiliation 
in a publication then we put a self-link with weight 1. 

In the citation network between countries (cities) 
nodes are papers which are linked if one paper cites the 
other. If a paper written by authors with n affiliations 
cites a paper written by authors with m affiliations we 
put nx m directed connections from each of the n citing 
countries (cities) to each of the m cited countries (cities), 
every link having weight l/(nm). The total weight of a 
directed link between two countries (cities) is the sum of 
all the weights over all the citations in the dataset. Since 
there can be multiple affiliations from the same country 
(city) in a publication, there are self-loops both in the 
world citation and in the world collaboration networks. 
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Appendix A: Methods 
1. Data description 



be used in our study. Further although the affiliations are 
being recorded with increasing consistency, their use still 
poses major challenges in uniquely and accurately iden- 
tifying them. For this reason, we parsed the affiliations 
of all publications and have determined the geographic 
location only at the city and country level. We also we 
use the publicly available resources (www. wikipedia. org 
and maps . google . com) to disambiguate the names of the 
places in case there are multiple name variation, typos 
and name changes during the time period of study. 



2. GDP 

The gross domestic product (GDP) is the value of all 
final goods and services produced within a nation in a 
given year and is the primary indicators used to gauge the 
health and size of a country's economy. We consider the 
average GDP (in US dollars) of a country during 2003- 
2010. A nation's GDP at purchasing power parity (PPP) 
exchange rates is the sum value of all goods and services 
produced in the country valued at prices prevailing in the 
United States. This is the measure most economists pre- 
fer when looking at per-capita welfare and when compar- 
ing living conditions or use of resources across countries. 



3. R&D spending 

Expenditures for research and development are cur- 
rent and capital expenditures (both public and private) 
on creative work undertaken systematically to increase 
knowledge, including knowledge of humanity, culture, 
and society, and the use of knowledge for new applica- 
tions. R&D covers basic research, applied research, and 
experimental development. 



For our study we used all publications in English in 
the databases of Science Citation Index Expanded, So- 
cial Sciences Citation Index, and Arts & Humanities Ci- 
tation Index for the years 2003-2010. The database of the 
Institute for Scientific Information (ISI) Web of Science 
also includes publications in other major languages, but 
consists of a relatively small number of items, accounting 
for < 5% of total publications. For each publication in 
the database, we have the name of the journal in which it 
is published, the volume and page number of the publi- 
cation, its year of publication, the names of the authors, 
the list of their affiliations and its references and other 
additional information. We used the list of references to 
construct the network of citations between papers. For 
each publication we extracted the city and country of 
authors institutions from the affiliation data. Whenever 
a publication has several authors, it is counted and as- 
signed to each location. Note that we only have the list of 
authors and the list of affiliations for each paper, however 
there is no corresponding match between these two lists 
and hence the individual level author affiliation can not 



4. Number of researcher 

Researchers in R&D are professionals engaged in the 
conception or creation of new knowledge, products, pro- 
cesses, methods, or systems and in the management of 
the projects concerned. Postgraduate PhD students en- 
gaged in R&D are included. 



5. Statistics 

To fit the data and calculate different estimates we use 
the following methods: 



a. Estimation of standard errors 

Bootstrapping is a distribution-free re-sampling 
method used to estimate the parameters of interest from 
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the empirical data. We have used this method in or- 
der to calculate the standard error of the mean. Let 
xi,X2, . . . ,x n be the dataset with mean x. The stan- 
dard error is then calculated as follows [41]: (i) Draw 
N samples each of size n with replacement from the 
original data, (ii) For each of the N samples calculate 
the sample mean a?i, xjv (hi) The standard error is 

then given by, SEE(x) = yj Y.iLi(%i ~ ^i) 2 > where 

Xi = TV -1 Xi is the mean of the TV bootstrap sample. 
In this study we have used 10 4 bootstrapped samples, 
i.e., N = 10 4 . 



b. Estimation of significance difference 

The above bootstrapping procedure however does not 
tell whether the difference in the means of two distribu- 
tions is significant or not. In this case the re-sampling 
has to be performed according to an appropriate null 
hypothesis, whereas for standard errors the re-sampling 
procedure was unrestricted. 

Let us consider two independent samples 
and yi, . . . , y m: and suppose that we are interested in the 
difference in the population means, S = x — y. Consider 
that the null hypothesis is Hq : x — y = 0. We create 
the bootstrap sample by choosing n elements without 
replacement from the pooled set x±, . . . , x n , y±, . . . , ?/ m . 
The remaining m elements constitute the other sample. 
We then calculate the mean of both these samples and 
determine the difference between them, say Si = Xi — yi. 
In analogous fashion N re-samples are made, and the 
bootstrap p value is defined as p = (#(^^^Vz))+i ^ j n 
this study we have used 10 4 — 1 bootstrapped samples. 



c. Power-law exponent 

We use maximum likelihood techniques to estimate the 
scaling exponent of power law distributions [42] . 



d. Regression Coefficient 

We used the linear regression analysis to study the re- 
lationship between the corresponding variables. We de- 
termine the regression coefficient using the ordinary least 
squares. The error term of the regression coefficient rep- 
resents the standard error of the estimate. 



6. Map construction 

Statistical data with embedded geographical informa- 
tion can be visualized with standard maps which are 
color coded by region. However these maps are some- 
times hard to interpret as the statistical measures are 



often correlated with the other indicators. We have 
used a diffusion-based method to create different density- 
equalizing maps [28]. In this method we start with an 
inhomogeneous distribution of the research contribution 
(in terms of citations, say) and let the diffusion process 
evolve until a homogeneous equilibrium state is reached: 
the displacements are then reinterpreted to generate the 
cartogram. 



Appendix B: Results 

We consider the research contribution of each coun- 
try in terms of the number of publications iVp u b, nor- 
malized by the number of participating countries in that 
publication. To visualize the results, we create a car- 
togram in which the geographic regions are deformed and 
rescaled in proportion to their relative research contribu- 
tion [28]. We observed that the contribution of different 
countries in terms of publications is heterogeneous and 
varies over 6 order of magnitude. Fig. |B.lft . shows that 
North America (32.4%), Europe(33.7%) and Asia(27.4%) 
have prominent contribution in terms of the number of 
publications. On the other hand, Africa, South America 
and Oceania contribute less than 7% of world's publi- 



cations. Table. B.l shows the contribution, number of 
countries and cities in each continents. It also indicates 
the statistics of the top countries of each continent. It is 
evident that the United States are the leading country in 
the world both in terms of publications and citations to 
them. It is followed by China, United Kingdom, Japan, 
and Germany in terms of publications, whereas in terms 
of citations it is followed by United Kingdom, Germany, 
Japan, and China. We indicate the fraction of total pub- 
lications /pub? the fraction of total citations received /cite 
and the average number of citations per paper, for coun- 
tries that received more than 0.005% of world citations. 
Countries are listed in decreasing order of the fraction 
of total citations received. The superscripts in /p u b and 
/cite indicate the world ranking of that country accord- 
ing to the numbers of publications and citations, respec- 
tively. We then consider the contribution in terms of the 
number of publications at the level of cities. In Fig. |B.lp 
we plot the probability distribution of the cities' contri- 
butions in terms of their publications and observed that 
it follows a power law scaling behavior with exponent 
1.45 ± 0.01. By plotting the out-degree against the out- 
strength, we find that there is po wer la w scaling behavior 
with (s out )(k out ) oc (k out ) a 



(Fig. 



B.l 



)). However, there 



are two distinct scaling regimes: for nodes with small kf 1 
(< 200) the exponent is a = 0.82±04, while for large k° ut 
(> 200) the exponent is a = 2.26±0.07. The super-linear 
behavior suggests that stronger links are more frequently 
connected to high out-degree nodes. 

Next we consider the average number of citations 
per paper of ea ch country and plot it on a colorpleth 
map (Fig. B.2A.). For calculating the average citation 
of a country we consider all its publications and count 
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Figure B.l. Research contribution in terms of number of publications. (A) Map of the country's research contribution, where 
the area of each country is scaled and deformed according to its number of publications. (B) The probability distribution 
function of the research contribution of cities in terms of their number of publications. The dashed line shows a power law 
scaling behavior with exponent 1.45 =b 0.01. (C) Node out-strength against its out-degree for city citation network. There are 
two distinct power law scaling regions, with scaling exponents 0.82 ± 0.04 and 2.26 =b 0.07 for low and high degree (> 200) 
nodes, respectively. 




Figure B.2. (A) Average citation of each country. World map where countries are color coded based on the average number 
of citations per publication. Most countries stay below the world's average of 7.67. (B) Weight of the links against the product 
of the strengths of the endpoints in the collaboration network of cities. There are two different scaling regions (2 x 10 7 ), with 
exponents 0.16 ± 0.01 and 0.92 ± 0.03. 



the total number of citations to all these articles during 
the period of 2003-2010. In the case where a publica- 
tion has multiple affiliations from different countries, it 
is counted multiple times for the countries' averages, once 
for each of the affiliated countries. In Table IB. 11 we have 



also given the average number of citations per paper of 
the top countries in each continent. The world average 
is 7.67. United States, Canada, Australia and most of 
the European countries have average number of citations 
larger than the world average. In Europe Switzerland 
leads the table, followed by Denmark and Netherlands. 
In contrast most of the countries from Asia stay below 
the world average, the only exception being Israel. Most 
of the countries in Africa and South America are below 
the world average as well. Other notable countries are 
Bermuda (16.97±5.95), Gambia (16.17±3.10), Panama 
(12.41± 0.68), Iceland (11.43± 0.71), Seychelles (ll.lli 
2.40), Guinea-Bissau (10.10± 0.97), Costa Rica (9.82± 
0.93), and Austria (9.75± 0.09). For the collaboration 



network of cities we plot the weight of the links against 
the product of the strengths of the connecting nodes, ex- 
pressing the expected weight of random collaborations 



(Fig. B.2^). As for citations we find that w%° 1 oc (siSj) a , 



with two different scaling exponents. If SiSj < 2 x 10 7 , 
a = 0.16 ± 0.01 (R = 0.11), whereas if SiSj > 2 x 10 7 
a = 0.92 ±0.03 (#=1.18). 

In Fig. |B.3| A,B we plot the probability of existence of 
a link as a function of the product of strength of the 
end-points of the link. We found that as the product 
increases, both in the collaboration and the citation net- 
work the probability of link existence increases, as ex- 
pected. In Fig. |B.3P ,D we show the variation of the 
link weight against the distance between the end-points. 
We found that both in the collaboration and the cita- 
tion network on the average the link weight decreases as 
a power-law with exponent 0.31 ± 0.01 and 0.22 ± 0.01, 
respectively. In this figure, while calculating the aver- 
ages we have only considered the existing links between 
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Figure B.3. Gravity law in the world collaboration and cita- 
tion networks. (A) Variation of the probability of existence of 
a link between two nodes as a function of the product of their 
strengths in the (A) collaboration network and (B) citation 
network of cities. Variation of the average link weight against 
the distance between the cities in the (C) collaboration net- 
work and (D) citation network. For each distance the average 
ratio is also shown. In this case only the existing links are 
considered while calculating the averages. The solid line indi- 
cates a power law behavior with exponent a = 0.31 ±0.01 and 
0.22 ± 0.01 respectively. Variation of the average link weight 
against the distance between the cities in the (E) collabora- 
tion network and (F) citation network. For each distance the 
average ratio is also shown. In this case all possible node pairs 
are considered in order to calculate the average, i.e., links that 
do not exist are considered with weight 0. The solid line in- 
dicates a power law behavior with exponent a = 0.88 =b 0.02 
and 0.51 =b 0.02, respectively. 



nodes. However, in the main text we have seen that the 
probability of link existence also decreases with distance. 
If we take this information while calculating the aver- 
ages, i.e., we consider the non-existent links by assigning 
weight zero to them, we found that in both the collabo- 
ration and the citation network, the average link weight 
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Figure B.4. Relation between research contribution in terms 
of number of publications and funding. Country's number of 
publications against the (A) expenditure in research and de- 
velopment (in million dollars, and purchasing power parity), 
(B) number of researchers in that country. The solid line in- 
dicates a scaling with exponent 0.90 =b 0.03 and 0.96 =b 0.03, 
respectively. (C) The plot of average spending per researcher 
against the average number of citation per paper of that coun- 
try. The average number of citations is now defined as the 
ratio of the normalized number of citations and normalized 
number of publications (see text). The horizontal line indi- 
cates world average, the vertical line indicates the spending 
of 120 000 $ per researcher. 



decreases with distance as a power law, with expo nent 
0.88 ± 0.01 and 0.51 ± 0.01, respectively (Fig. |R3^,F). 
Note that this property is different from what has been 
observed in the mobile phone communication network, 
where it was shown that the weight of the existing links 
are independent of the distance, whereas the overall link 
weight decrease as a result of decreasing probability of 
having a link as the distance increases [2]. 

In the main paper we have considered the research per- 
formance of each country based on the number of cita- 
tions. In addition, here we consider the performance of 
a country based on its number of publications. As be- 
fore, in Fig. |B.4| A, we plot the research contribution in 
terms of the number of publications iVp u b against the 
countries' R&D expenditure in terms of purchasing power 
parity (PPP). We found that this indicator also scale al- 
most linearly with the spending. We next consider the 
dependence of research performance on the number of re- 
searchers in that country (Fig. 



B.4B). The research con- 



tribution in terms of publications also scale linearly with 
the number of researchers in that county. 

Finally as a measure of the average publication qual- 
ity of a country we consider the ratio of the normalized 
number of citations and normalized number of publica- 
tions of that country. This is an alternative measure 
of the average number of citations per paper we men- 
tioned above, which is not normalized by the number of 
authors in a paper. In the previous measure each pub- 
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lication from a country (independently of the number of 
participating countries) gets equal weight while calculat- 
ing the average. In this other measure, if there are n 
countries in a publication, each country would get 1/n 
as credit for that publication, so that publication would 
give a lower contribution to the average number of cites 
per paper than before. In Fig. |B.4p we plot the new 
quantity against the average spending per researcher of 



the country (R&D expenditure divided by the number of 
researchers). Although this plot is similar to the one in 
the main paper, there are certain differences in the aver- 
age number of citations of some countries. For example, 
Italy, Spain, Norway are now below the world average. 
This means that the publications from these countries 
with international collaborators contribute significantly 
to the average impact of their scientific production. 



Table B.l: Research contribution of different continents and their top countries. The number of countries and cities in each 
continent are indicated by A/countries and iVcities, respectively. Fraction of publications /p u b, fraction of citations received /cite 
and the average number of citations per paper of each continent is also indicated. For top countries in each continent we list 
the fraction of publications /p u b, fraction of citations received /cite, the average number of citations per paper. The superscript 
indicates the countries' rank in the world in terms of number of publications and number of citations. Only countries that 
receive more than 0.005% of all citations are shown. 



Continent ^Countries ^Cities /pub /cite Avg. Cites Country / Pub /cite Avg. Cites 

(in %) (in %) name (in %) (in %) 



Africa 



57 



749 1.32 0.65 5.00±0.05 



South Africa 
Egypt 
Tunisia 
Nigeria 
Kenya 
Morocco 
Algeria 
Tanzania 
Uganda 
Cameroon 
Ethiopia 
Ghana 
Zimbabwe 
Malawi 
Senegal 
Botswana 
Gambia 
Cote d'lvoire 



0.430 33 
0.286 38 
0.100 52 
0.126 5 " 
0.038' 



65 
-60 



0.248 37 
0.128 40 
0.036 52 
0.031 56 
0.028 58 



0.055 DU 0.025' 
0.067 54 
0.019 83 
0.018 85 
0.021 77 
0.022 75 
0.016 86 
0.011 95 
0.009 103 
0.008 104 
0.010 97 
0.003 128 
0.006 107 



60 

0.023 62 
0.014 74 
0.014 75 
0.010 80 
0.009 81 
0.008 85 
0.007 87 
0.006 88 
0.006 90 
0.005 95 
0.005 96 
0.005 97 



5.92±0.08 
3.78±0.05 
3.33±0.13 
2.82±0.25 
7.55±0.29 
4.20±0.10 
3.01±0.08 
7.27±0.29 
7.04±0.27 
4.72±0.18 
4.36±0.17 
5.29±0.21 
5.92±0.28 
7.11±0.32 
6.72±0.29 
5.46±0.54 
15.87±3.02 
7.20±0.41 



Asia 



49 



3853 27.36 17.71 5.58±0.01 



Japan 


6 


,457 4 


5 


.939 4 


7. 


.68±0. 


,03 


China 


7 


.216 2 


4 


.304 5 


5. 


,05±0. 


.02 


South Korea 


2 


.509 10 


1 


.582 13 


5. 


.38±0. 


.04 


India 


2 


:.727 9 


1 


.398 15 


4. 


35±0. 


,03 


Taiwan 


1 


.671 15 


1 


.037 16 


5. 


19±0. 


.04 


Israel 





.863 22 





.837 20 


8. 


86±0. 


,10 


Turkey 


1 


.450 17 





.667 22 


3. 


.89±0. 


,03 


Russia 


1 


.875 13 





.622 24 


3. 


92±0. 


,05 


Singapore 





.521 29 





.461 28 


7. 


,29±0. 


,08 


Iran 





747 23 





.308 31 


3. 


31±0. 


,03 


Thailand 





.244 41 





14? 38 


5. 


.88±0. 


,17 


Malaysia 





.195 42 





.069 45 


3. 


.22±0. 


,07 


Pakistan 





.175 44 





.059 48 


3. 


23±0. 


,07 


Saudi Arabia 





.131 49 





.046 50 


3. 


12±0. 


,07 


Jordan 





.061 58 





.023 61 


3. 


.27±0. 


,11 


Vietnam 





.037 68 





.020 66 


5. 


59±0. 


,30 


Indonesia 





.032 70 





.019 67 


5. 


73±0. 


,23 


Kuwait 





.044 61 





.018 68 


3. 


,80±0. 


,16 


Bangladesh 





.041 63 





.018 69 


4. 


74±0. 


,17 


Lebanon 





.038 66 





.018 70 


4. 


,47±0. 


.13 


UAE 





.041 62 





.018 71 


4. 


03±0. 


,14 


Philippines 





.035 69 





.017 72 


6. 


20±0. 


,30 


Cyprus 





.027 73 





.012 76 


4. 


,31±0. 


,15 


Sri Lanka 





.022 76 





.Oil 79 


4. 


,90±0. 


,18 


Armenia 





.026 74 





.009 82 


6. 


18±0. 


,31 


Oman 





.021 79 





.008 86 


3. 


,50±0. 


,14 


Georgia 





.020 82 





.006 89 


3. 


,94±0. 


.20 


Nepal 





.012 92 





.006 91 


5. 


36±0. 


,29 


Uzbekistan 





.020 81 





.006 92 


3. 


,50±0. 


,17 
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Table B.l - Continued from previous page 



Continent ^Countries Amities huh /cite Avg. Cites Country / Pub 

(in %) (in %) name (in %) 



fate Avg. Cites 
(in %) 



Europe 



47 



6625 33.69 35.25 9.29±0.01 



United Kingdom 
Germany 
France 
Italy 
Netherlands 

Spain 
Switzerland 
Sweden 
Belgium 
Denmark 
Finland 
Austria 
Poland 
Norway 
Greece 
Portugal 
Czech Republic 
Ireland 
Hungary 
Slovenia 
Ukraine 
Romania 
Slovakia 
Croatia 
Serbia 
Bulgaria 
Estonia 
Lithuania 
Iceland 
Belarus 
Latvia 
Luxembourg 
Moldova 



6.509 3 
5.131 5 
3.611 7 
3.415 8 
1.829 14 
2.482 11 
1.114 19 
1.227 18 
0.923 21 
0.655 25 
0.640 26 
0.595 27 
1.110 20 
0.506 30 
0.684 24 
0.460 32 
0.464 31 
0.340 36 
0.338 37 
0.171 45 
0.271 40 
0.274 39 
0.150 48 
0.164 47 
0.167 46 
0.125 51 
0.063 57 
0.095 53 
0.030 71 
0.064 56 
0.021 78 
0.012 91 
0.011 96 



7.453 2 
6.299 3 
4.034 6 
3.258 8 
2.331 9 
2.258 11 
1.600 12 
1.436 14 
1.004 17 
0.838 19 
0.672 21 
0.653 23 
0.579 25 
0.483 26 
0.468 27 
0.345 30 
0.301 32 
0.298 33 
0.251 36 
0.096 41 
0.088 42 
0.079 43 
0.072 44 
0.068 46 
0.061 47 
0.056 49 
0.041 51 
0.035 53 
0.031 57 
0.020 65 
0.009 83 
0.008 84 
0.006 93 



9.91±0.04 
10.41±0.04 
9.67±0.04 
8.59±0.04 
11.08±0.07 
8.09±0.05 
12.38±0.09 
10.59±0.09 
10.02±0.08 
11.45±0.12 
9.59±0.10 
9.75±0.11 
5.43±0.05 
8.98±0.10 
6.30±0.06 
6.93±0.08 
6.31±0.08 
8.20±0.16 
7.61±0.13 
5.41±0.09 
3.46±0.07 
3.30±0.09 
5.12±0.12 
4.53±0.12 
3.42±0.07 
4.76±0.09 
6.91±0.21 
3.91±0.13 
11.46±0.63 
3.37±0.10 
5.34±0.43 
6.58±0.34 
4.94±0.31 



North America 37 



5346 32.40 42.33 10.36±0.02 



United States 
Canada 
Mexico 
Puerto Rico 

Cuba 
Costa Rica 
Panama 



28.116 1 
3.616 6 
0.523 28 
0.037 67 
0.040 64 
0.014 88 
0.009 102 



38.216 1 
3.728 7 
0.292 34 
0.028 59 
0.022 63 
0.012 77 
0.012 78 



10.67±0.02 
9.15±0.05 
5.57±0.10 
7.66±0.26 
4.81±0.14 
9.93±0.87 
12.43±0.81 



Oceania 



21 



844 2.89 2.67 8.22±0.05 



Australia 
New Zealand 



2.448 12 2.301 10 
0.425 34 0.354 29 



8.36±0.05 
7.60±0.10 



South America 14 



782 2.34 1.39 5.75±0.04 



Brazil 
Argentina 

Chile 
Colombia 
Venezuela 
Uruguay 

Peru 
Ecuador 



1.551 lb 
0.399 35 
0.193 43 
0.066 55 
0.060 59 
0.027 72 
0.019 84 
0.008 105 



0.871 18 
0.261 35 
0.136 39 
0.034 54 
0.034 55 
0.021 64 
0.014 73 
0.005 94 



5.21±0.04 
6.31±0.10 
7.42±0.16 
5.65±0.19 
6.11±0.28 
6.81±0.21 
7.68±0.30 
7.39±0.37 



